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Abstract 

We study the elliptic flow V2 of multi-strange hadrons such as the (p, H and as 
a function of transverse momentum in the recombination and fragmentation model 
and compare to a standard hydrodynamic calculation. We find that the measure- 
ment of V2 for the (f) and will allow for the unambiguous distinction between 
parton recombination and statistical hadro-chemistry to be the dominant process 
in hadronization at intermediate transverse momenta. 
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1 Introduction 



Elliptic flow of liadrons - more precisely defined as the second Fourier co- 
efficient of their azimuthal momentum distribution - has been suggested as 
a sensitive probe for the buildup of pressure in the early reaction stage of 
relativistic heavy- ion collisions [1]. Data from Au-|-Au collisions at the Rel- 
ativistic Heavy Ion Collider (RHIC) has revealed unexpected features in the 
behavior of the elliptic flow V2 as a function of transverse momentum Pt'- for 
non-central collisions V2 first rises and then saturates at higher Pt [2,3,4,5,6], 
the saturation occurring at higher Pt and at higher values for V2 for baryons 
than for mesons. 
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It was first pointed out by Voloshin that this surprising behavior can be ex- 
plained in a naive parton recombination model for hadronization [7]. Several 
more detailed studies on elliptic flow based on recombination have been carried 
out subsequently [8,9,10,11,12,13]. In addition to elliptic flow, parton recombi- 
nation has been very successful in explaining the behavior of particle spectra 
and ratios at intermediate transverse momenta. The general theory of recom- 
bination as a hadronization mechanism in dense parton systems is described 
in great detail in [8,9,13,14,15,16,17,18,19]. In this letter we present quantita- 
tive results for the elliptic flow of multi-strange particles in the recombination 
plus fragmentation approach - the theoretical foundation being based on our 
previous work in [9]. 

At low Pt the experimental data show that the elliptic flow of a heavy particle 
is smaller than that of a light particle. This mass-dependence can be under- 
stood in hydrodynamic models where particles with higher mass experience 
the collective movement of the medium less than lighter species [20,21]. How- 
ever, the saturation of V2 above Pt ~ 1.5 GeV/c and its dependence on the 
hadron species in this region contradicts a purely hydrodynamic model. 

At high Pt, the energy loss of fast partons in the dense medium has been pro- 
posed as a mechanism that can translate the anisotropy in the initial Au+Au 
system into an momentum anisotropy in the flnal state as well [22,23,24]. The 
observation of jet quenching at RHIC energies is a strong indication for the ex- 
istence of partonic energy loss which may lead to such an azimuthal anisotropy 
[25]. 

In the most interesting region of intermediate Pt, between 2 and 5 GeV/c, 
the use of hydrodynamics becomes questionable, since the fundamental as- 
sumption of an inflnitely small mean free path is not valid anymore. On the 
other hand, the use of perturbative QCD is not well justifled either, since the 
relevant scale ~ Pt is rather small and may lead to uncontrollable higher twist 
corrections. In addition, the strong energy loss further suppresses the contri- 
bution from perturbative parton fragmentation. These considerations suggest 
recombination as the dominant hadronization mechanism in the intermediate 
Pt region. 

To clarify the situation and to establish the dominant mechanism for hadroniza- 
tion and the creation of hadron V2 in the intermediate Pt domain we propose 
the measurement of w 2 for multi-strange hadrons, in particular for the (p and 
the Q. In a hydrodynamic picture the differences between hadron species are 
predominantly driven by the mass differences. In the recombination formalism 
the mass only plays a role at low Pt, while strong deviations appear at inter- 
mediate Pt due to the different valence structure of baryons and mesons. Since 
the mass of the 4> is very close to that of the A, a statistical hadronization 
scheme as employed by hydrodynamics would predict the (j) to behave similar 



2 



to the A, while recombination imphes that it behaves similar to other mesons, 
e.g. the kaon. 



2 Recombination plus fragmentation approach 

In [8,9] we found that hadron production at intermediate Pt is governed by an 
interplay of fragmentation and recombination. Recombination is more efficient 
than fragmentation if the parton spectrum is exponential, while fragmentation 
dominates for the (perturbative) power law tail of the partonic Pt spectrum. 
E.g. for mesons M recombining on a hypersurface E, the momentum distri- 
bution is given by 

E^ = Cm j^daa^^^^ jyxWa{R; xP+)\<Pm{x)\'w,{R; {l-x)P+), (1) 

where 4>m{x) is a wave function for the meson on the light cone, x is the 
momentum fraction of valence quark a, the w are the momentum distributions 
of the quarks a and b and Cm is a degeneracy factor. We refer the reader to 
[9] for more details. 

The parton momentum distribution at hadronization is an input in Eq. (1). 
For quantitative predictions we assume a thermal distribution, 

where v{R) is a four velocity and 7^ is fugacity factor for each parton species 
a. The longitudinal and transverse spatial distributions are described by the 
width A for the space-time rapidity rj and a profile function /(p, 0) for trans- 
verse size and azimuthal angle p and (f). We set A = 2. f{p,(l)) is impact 
parameter dependent and is scaled with the volume of the overlap zone of two 
nuclei in non-central collisions. 

We assume that the hadronization process takes place at a temperature T = 
175 MeV on a hypersurface S given by r = \Jt'^ — = 5 fm/c. The radial 
flow velocity vt is fixed at 0.55c for all values of p and for all impact param- 
eters. These input parameters were determined in [9] and are consistent with 
the observed hadron spectra and ratios from PHENIX and STAR. The re- 
combination calculation has to be supplemented by a pQCD calculation using 
fragmentation functions. However, no information is available concerning the 
fragmentation functions of strange hadrons beyond kaons and As. 
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3 Elliptic flow 



The azimuthal asymmetry V2 is defined by 

where $ is the azimuthal angle in momentum space. 

At low and intermediate Pt the initial geometric anisotropy results in an 
anisotropic pressure gradient that implies an elliptic flow profile. This is para- 
metrized via an anisotropic parton momentum distribution at the time of 
hadronization at which partons then hadronize via recombination. We use Eq. 
(2) together with an elliptic modulation of the transverse rapidity 

77T(0,Pr) =r^^(l -i^(pT)cos20) (4) 

to determine V2 for the partons at hadronization. 77^ is the rapidity given by 
the radial flow velocity vt — 0.55c, so that tanhr^y = 0.55. This is inspired by 
a hydrodynamical description of the parton phase at low pt- Hydrodynamics 
works well to describe the measured V2 for hadrons up to Pt = 1.5 GeV/c 
[4,26,27,28]. In order to accommodate the deviation from ideal hydrodynamics 
at higher Pt we assume that F{pt) takes the from 

F{Pt) = .2 ^ (5) 

1 + [pt/PoY 

where for a given impact parameter h 

w{h) + l(h) ^' 

is determined by the collision geometry with width w{h)—RA — b and length 
l{b) = ^JR\ - (6/2)2 of the collision zone. Ra is the radius of a gold nucleus. 

We assume that V2 = V2 = V2 = V2 and f| = f|. The slight difference between 
light quarks and strange quarks originates from the mass difference. We use 
constituent masses mu,d = 260 MeV and = 460 MeV. The parameter 
Po — 1.1 GeV/c was obtained in [9] by a fit to the eUiptic fiow of pions 
obtained by PHENIX [2] and is also consistent with data on V2 for charged 
hadrons, protons, kaons and As [9]. 

In [9] we showed that for recombining mesons M and baryons B the anisotropy 
V2 can be written as 

M/p X ^ J dx\<PM{x)\^[v2ixPT) + t'2((l - x)PT)]kM{x, Pt) . . 

^ Jdx\Mx)\'[l + 2vl{xPTH{{l-x)PT)]kM{x,PTy ^ ^ 
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where kui^, Pt) and Pt) are given by 



kM{x, Pt) = Ki 



{^Jml + x^Pl + ^ml + (1 - x)2P, 



, (8) 



(9) 



cosh r^T- 



E 

j=a,b,c 



ml + xlP^ 



(10) 



a, h and c stand for the valence partons respectively. We use the short notation 
/ Vxi = Jq dxadxf,dx,,5(xa + Xh + Xc — 1) in Eq. (8). Using infinitely narrow 
delta shaped wave functions |0p ~ 6 that distribute the momentum of the 
hadron equally to all valence quarks and neglecting higher order terms one 
obtains the simple scaling law [5,6,7,9,11,12], 



v^(Pt) = nv^, ( -Pt 



where n is the number of valence quarks in the hadron. 



(11) 



At high Pt the elliptic fiow is dominated by the contribution of fragmentation 
and its value is determined by the azimuthal anisotropy due to jet energy 
loss. Since fragmentation functions of multi-strange particles are basically un- 
known, a precise calculation is not possible. However, the energy loss and hence 
the resulting anisotropy are partonic in nature and one can argue that in the 
ratio of Eq. (3) the effects of fragmentation functions tend to cancel. Indeed, 
we have shown in [9] that for the part of V2 coming from fragmentation the 
difference between different hadron species is numerically negligible. Keeping 
this in mind we approximate the component of f 2 coming from fragmentation 
for baryons and mesons by that of As and kaons respectively. 



The total eUiptic fiow can be written as 

V2{Pt) = r{PT)v2,R{PT) + (1 - r(PT)KF(PT), 



(12) 



where r{PT) is the ratio recombination / (recombination + fragmentation) for 
the hadron spectrum as a function of transverse momentum. 



4 Numerical results 



In the following we shall discuss our results for the elliptic flow of (multi- 
) strange particles. All results shown are for Au-|-Au coUisions at ^/s — 200 GeV 
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Fig. 1. Elliptic flow of u quarks and s quarks as a function of pT (thermal phase 
only). 

with an impact parameter b of 8.0 fm. Starting point of our calculation is a 
system of constituent quarks prior to hadronization. Figure 1 shows f 2 of u and 
s quarks as a function of pt from Eqs. (3,4,5). The slight difference between 
u and s quarks stems from their different masses. The mass effect disappears 
at high Pt- The assumption of equal elliptic flow of s and u, d quarks is well 
justified given the nearly identical behavior of the kaon and pion elliptic flow 
as a function of pt, as can be seen in [9]. Note that a different value for the s 
quark elliptic flow will lead to a characteristic flavor ordering for the elliptic 
flow among multi-strange mesons and baryons respectively, depending on their 
strange quark content [10,11,12]. Such an effect has not been observed in the 
data so far. 



Figure 2 shows V2 for multi-strange mesons (top frame) and baryons (bottom 
frame) as a function of Pt- The experimental data are taken from the STAR 
collaboration [6,29]. Thick lines denote the full recombination-|-fragmentation 
calculation (labeled as R-fF in the figure), whereas thin lines refer to recom- 
bination from thermal quarks alone (labeled as R in the figure). Under the 
assumption that the s quark elliptic flow is identical to u, d quark elliptic flow 
all mesons show the same behavior. The V2 for is identical to that of the K- 
This is most prominent in the saturation region of V2- In the recombination 
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approach the value of f 2 there is solely determined by the number of valence 
quarks of the hadron and the partonic V2- The mass effect which drives the V2 
characteristics in hydrodynamical calculations is negligible for the saturation 
feature. The same systematics can be observed for baryons as well, leading to 
nearly identical elliptic flow values for A, S and Experimental data conflrm 
our calculation in the R+F approach for kaons, As and S between 1.5 and 6 
GeV/c. Wc have shown in [9] that the dependence of elliptic flow on the shape 
of the wave function used in Eqs. (7,8) is negligible for transverse momenta 
larger than 2 GeV. For further studies on the influence of the wave function 
see e.g. [15]. 

At Pt lower than 1.5 GeV the data exhibit the well known splitting of the 
V2 curves due to the mass of the hadrons. This is only poorly reproduced by 
the recombination formalism since binding energies are not properly taken 
into account, whereas hydrodynamic calculations perform very well in this 
domain. 

The fragmentation contribution to V2 and the ratio t^Pt) for kaons and As 
are taken from [9] . As discussed above we use the values for kaons to estimate 
both V2 from fragmentation and r{PT) for the 0. Moreover we use the values 
for the A to estimate the corresponding quantities for the S and the Q in the 
fragmentation process. The tendency of fragmentation functions to cancel in 
Eq. (3) makes this a good approximation and also makes the difference between 
mesons and baryons vanish. This implies a universal flat behavior for > 6 
GeV/c where recombination is suppressed and the azimuthal anisotropy is 
dominated by pQCD processes and parton fragmentation. The merging of 
meson and baryon V2 into one universal curve at high Pt can be seen even 
better in Fig. 3. Multistrangc hadrons with high masses might be disfavored 
by the fragmentation process. We expect the relative weight of recombination 
to be even larger for multi-strange particles: r^{PT) > rxiPr), tq.{Pt) > 
fAiPr) > fsiPr)- Therefore hadrons like the S and the fl might reach the 
universal curve for V2 from fragmentation even later than indicated in Fig. 2 

Figure 3 elucidates the systematic differences between V2 originating from sta- 
tistical hadronization in a hydrodynamic scenario (blast wave model [21,30]) 
vs. hadron V2 generated from a recombination-|-fragmentation scenario: thin 
lines represent a hydrodynamic calculation of the elliptic flow of 0. A and Q 
as a function of Pt- Since the masses of the and the A are similar, their f 2 is 
identical as well in the hydrodynamic calculation. Due to its much larger mass 
the Q exhibits a smaller V2, clearly distinct from that of the other two hadron 
species. The situation reverses in the region where recombination dominates: 
there the number of valence quarks in the hadron determines the resulting V2, 
leading to a distinct difference in the f 2 between the and the A, which now 
coincides with the Q. At highest Pt all curves merge to the universal V2 from 
fragmentation. 
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Fig. 2. Elliptic flow for multi-strange mesons (top frame) and baryons (bottom 
frame) in the recombination-l- fragmentation approach compared to data from STAR 
(symbols) (A + A [6], H + H [29]). Thick lines denote the full calculation, whereas 
thin lines refer to the recombination contribution. 
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Fig. 3. Elliptic flow for some strange hadrons in the recombination+fragmentation 
approach (thick lines) and in a hydrodynamic calculation (thin lines). 

We assume a sudden freeze-out and neglect interactions in the hadronic phase 
in both descriptions, recombination plus fragmentation and hydrodynamics. 
This seems to be reasonable for particles with high and intermediate transverse 
momentum and is supported by the success of the scaling law Eq. (11). The 
study of resonance particles is particularly important in this case [17,31]. We 
will investigate this in a forthcoming publication. 



5 Conclusion 



In brief, we investigated the elliptic flow of multi-strange hadrons in the re- 
combination plus fragmentation approach. We find that the behavior of f 2 as a 
function of Pt in this approach is dominated by the number of valence quarks 
of the respective hadron, leading to very similar values of the elliptic flow 
for all mesons and likewise for all baryons, nearly independent of the hadron 
mass, as opposed to conventional statistical hadronization in hydrodynamic 
calculations. In particular we find that f 2 for mesons is nearly identical to 
that of kaons above 2 GeV/c. This uses a parton V2 for strange quarks that 
is for all practical purposes equal to the V2 of the light flavors, in accordance 
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with RHIC data. The measurement of and fl eUiptic flow will thus permit 
to distinguish between statistical hadronization as employed by standard hy- 
drodynamic calculations vs. hadron production through parton recombination 
at intermediate Pt- At high Pt we flnd the elliptic flow to be dominated by 
fragmentation, leading to a universal curve above 6 GeV/c. 
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